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Abstract 



We present fermi liquid expressions for the low temperature behavior of 
the superfluid stiffness, explain why they differ from those suggested recently 
by Lee and Wen, and discuss their applicability to data on high-T c super- 
conductors. We find that a consistent description requires a strong, doping 
dependent anisotropy, which affects states near the zone corners much more 
strongly than those near the zone diagonals. 
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The work reported here was motivated by a recent paper of Lee and Wen on " Unusual 
Superconducting State of Underdoped Cuprates"0. Their point of view, which we share, is 
that because the superconducting state of the high-T c materials is apparently conventional 
and in particular supports well-defined quasiparticle excitations, it is appropriately described 
in terms of fermi liquid theory. Their main conclusion, that the temperature dependence 
of the superfluid stiffness is anomalous and implies unusual properties, is important and 
interesting. Their specific results must be interpreted with caution because their formulation 
of fermi liquid theory in the superconducting state (Eqs 1-4 of their paper) disagrees in an 
important way with the standard formulationi. In this note we present the usual fermi liquid 
theory of the charge response of the superconducting state, explain why it differs from Eqs 
1-4 of Lee and Wen, and apply it to penetration-depth data and discuss the underlying 
physics. 

We begin, however, with data. We write the penetration depth A, in terms of a superfluid 
stiffness per Cu02 plane p s via p s = h 2 c 2 b/167ie 2 \ 2 . Here b is the mean interplane spacing. 
This definition yields the coefficient of the |(W) 2 term in the Ginzburg-Landau theory 
but differs by a factor of 4 from the quantity defined in0. The most extensively studied 
family of materials is the 'YBCO' family; to avoid complications due to the C11O2 chains 
we consider the penetration depth for current flow in the basal plane but perpendicular 
to the chains. In YBCO the length b ~ 5.5A. The measured!! penetration depths are 
\ aa = 1600 for YBCI2CU3O7 and A aa = 2000 for YBCI2CU3OQS implying p saa = llmeV (O7) 
and p S aa = 7meV (Oqs)- In a clean weak-coupling superconductor p s may be determined from 
the conduction band plasma frequency u p via p s = (hu p ) 2 b / (16n e 2 ); the value calculated 
using band theory is@ (hu p ) 2 = 9eV 2 implying p sa a-band = 70meV. The renormalization of 
p s from the band theory prediction is large (factor of 6-10) and doping dependenti. 

The temperature dependence of p s for YBCO is given ini. For both O7 and Og.6 Ps is 
well fit by p s (T) = p s (T = 0)(1 — 0.5T/T C ). Converting the data to absolute units we 
find limj^o dp s /dT = 0.06meV/K for both materials. The theoretical result^ is dp s /dT = 
(ln2/2\/2)(vp/a)/dA/d9 where dA/d9 is the angular derivative of the gap function at the 
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gap node and a = 3.8A is the lattice constant. Photoemission results^! on BSCCO show 
that if the gap is defined as the displacement of the midpoint of the leading edge of the 
photoemission spectrum from the fermi level, then dA/d8 ~ AOmeV ~ A50K; using the 
band theory vp ~ 4eV — A yields dp s /dT = 0.5meV/ K. 

We now interpret these data in terms of fermi liquid theory. Fermi liquid theory is 
specified by the quasiparticle velocity vp{9) and the Landau interaction function f(9,9'); 
these are functions defined on the fermi surface which for the electronically two-dimensional 
high-T c materials is a simple closed curve Pf(9) with 9 the usual polar angle (we neglect the 
chain bands and any bonding-antibonding splitting arising from the bilayer structure). 

A basic assumption of the fermi liquid analysis is that the low-T dependence is due to 
thermally excited quasiparticles. Several authors have claimed^ that phase fluctuations 
of the superconducting order parameter can also lead to such a T-dependence. However, 
these models require a large density of normal electrons at T = to screen the Coulomb 
interaction. In high-T c materials there are no such electrons at T = 0; the phase fluctuations 
will therefore be gapped, giving an exponentially small contribution. Even if screening is 
included, we believe the models of refsiS yield a T 2 temperature dependence. Claims to the 
contrary seem to have arisen from Eq 7 of rei§ in which a term involving the Matsubara sum 
Y^n=°?oo was written as 2 * Y^^Lq so that the n = term was counted twice. Evaluation of 
Y^n=°?oo yields a T 2 temperature dependence 

The fermi velocity may be determined from the dispersion of quasiparticle peaks ob- 
served in photoemission. The interpretation of photoemission data on high-T c materials is 
still controversial, but there is general agreement that for momenta near the zone diagonal, 
well-defined peaks exist from which a dispersion can be extracted. The most detailed studies 
have focused on the BSCCO family of materials; the data for YBCO seem quite similar!. 
Fig. 3 of re shows the dispersion along the zone diagonal for overdoped, optimally doped, 
and underdoped BSCCO samples. The data imply a doping-independent quasiparticle ve- 
locity v = 1.3eV — A; this is reduced from the band-structure value {yp ~ 4eV) by about a 
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factor of 3, presumably because of electron-electron interaction effects. It is conventionally 
assumed (and has been demonstrated in some models0'lll) that the many-body renormal- 
izations which change the zone-diagonal velocity by a factor of three are relatively isotropic 
around the fermi surface, so we might expect that the velocity everywhere to be changed 
by about the same factor. This, combined with the observation that the charge response 
should be dominated by the zone diagonals where the velocity is largest, suggests that if 
the Landau parameter contribution were negligible p s would be renormalized by about a 
doping-independent factor of three from the band theory prediction. The smaller value and 
especially the doping dependence p s therefore suggests that the Landau interaction function 
plays an important role as predicted theoreticallySS. We now show that while the Landau 
function may be important, the conventional picture makes predictions for dp s /dT which 
are inconsistent with the data. 

The fermi liquid theory of the superconducting state was developed irM The result may 
be written: 

Psa b = ^ / de 1 d9 2 N{e 1 )N{e 2 )v a {e 1 )L{e 1 )T b {e 1 , e 2 ) (1) 

Here N(9) = \/v{9) is the quasiparticle density of states at position 9 on the fermi 
surface, v a is a component of the quasiparticle velocity and L{9) = nT J2 n A(9u> n ) 2 /[u^ + 
A(9, uj n ) 2 } 3 / 2 is the usual angle-dependent Yoshida function. The vertex T obeys the equation 

T = v + (fLT) (2) 

Here / is the Landau function defined above and products inside angle brackets are 
shorthand for integrals over the fermi line, weighted by density of states as in Eq. |l|. Eq. |2| 
may be solved formally and the solution used in Eq. |1|. Noting that at T = L = 1 and 
that the leading T-dependence is in L one finds: 

Psaa(T = 0)=<V a (l-f)- 1 V a > (3) 



lim dp saa /dT =< v a {l - f)- l dL/dT{l - f)- l v a > (4) 



Eq. [| differs from that written in refS by the two factors of (1 — / ) 1 . The authors of§ ob- 
tained their result by analogy to Leggett's theory of the spin susceptibility of 3 He; the results 
oi show the analogy is not appropriate. The physics is this: Landau parameters enter ex- 
pressions for physical quantities because a perturbing field H excites quasiparticle-quasihole 
pairs at all points on the fermi line; i.e. it modifies the underlying fermi distribution and 
because of the electron-electron interaction this changes the quasiparticle energy dispersion 
to order H, and thus affecting the number of excitations. This effect must be treated self- 
consistently. Now in 3 He the condensate does not respond at linear order to a magnetic field 
(at least in the low-T "BW" phase), and because of the superconducting gap the number of 
quasiparticles at low T is very small; thus at low T changes to the quasiparticle dispersion 
may be neglected and Landau parameters drop out. By contrast a superconducting con- 
densate does respond to a magnetic field (by producing a super cur rent); this supercurrent 
does correspond to a displacement of the underlying fermi distribution which does affect 
the superconducting quasiparticles, causing Landau parameter effects in the temperature 
dependence. 

In fermi liquid theory one defines the Landau parameters A\ a , F ls via < v a (l — 
f)~ l v a >=< v 2 a > /(l - A u /2) =< v 2 a > (1 + F ls /2). If the renormalization of the velocity 
is isotropic over the fermi surface we may evaluate < v 2 > by dividing the band theory p s by 
three (recall the density of states factor in the definition of <>) to obtain < v\ >= 23meV. 
By comparing this to the data we conclude that (l + F ls /2) ps 1/2 (0 7 ) and (l+F ls /2) 1/4 
(06.6)- Now in a Galilean-invariant fermi liquid Vp = pp/m* and < v\ >~ n/m* . In a two 
dimensional Galilean-invariant system the Landau parameter obeys 1 + F± s /2 = m*/m so 
the Landau parameter cancels the mass renormalization. In non-Galilean-invariant systems 
such as the high-T c superconductors there is no a-priori relation between F ls and m*/m. As 
we see, in the actual materials the observed p s is smaller than one would expect from the 
observed v ; thus Fi. s acts in the same direction as m*, instead of compensating it. 

In many fermi liquids the angular dependence of the interaction is not too strong, so 
one expands the operator f in appropriate harmonics and retains only the one most closely 
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corresponding to vv; this amounts to replacing the operator (1 — f)^ 1 by the scalar l + Fi s /2 
in both eqs. |3] and |J We refer to this as the "conventional fermi liquid" approximation. A 
particular realization of such a fermi liquid is found in the slave boson calculations o flfi To 
compute dp s /dT in this conventional fermi liquid approach we take the "band theory" value 
of 0.5meV/K obtained earlier, divide by a factor of three for the velocity renormalization 
and a further factor of 4 for Y Ba2Cu^Oy and 16 for Y Ba2Cu^O^,^i yielding dp s /dT = 
0. 04meV/ K for O7 and O.OlmeV/ K for 06.6- Although the estimate of dp s /dT for O7 is 
within a factor of two of the measured value of 0.06meV/K, the estimate for 06.6 is badly 
off. Similarly, reM found p s {T) =< vL(T)v > /(l + c < vL(T)v > with c w -A u ~ l/x, 
so dp s /dT ~ p s (T = 0) < vdL/dTv > jc ~ l/x 2 . Thus both the conventional fermi liquid 
approach and explicit calculations predict a strong doping dependence where none exists 
experimentally. In re^l the disagreement of the result of re fO with the observed dp s /dT 
was argued to indicate the inadequacy of the U(l) gauge theory formalism used in 0. The 
present analysis makes it clear that the difficulty is more general. 

The doubtful point of the discussion given above is the assumption of weak angular 
dependence around the fermi surface. This enters the argument in two places: first, that the 
factor-of-three velocity renormalization observed for quasiparticles along the zone diagonal 
applies all over the fermi surface, and second that the Landau interaction operator (1 — f)^ 1 
could be replaced by the Landau parameter 1 + Fi s /2. These two issues are clearly related: 
an interaction giving a relatively larger velocity renormalization at the zone corners than at 
the zone diagonals will lead to a Landau function with a stronger angular dependence. 

If the velocity renormalization increases rapidly as one moves away from the zone diagonal 
then the estimate of the renormalization of p s relative to band theory increases and the value 
of the Landau parameter decreases. There is evidence of anomalous flatness of bands and 
shortness of lifetimes in the vicinity of the zone corners^. In order to explain the data the 
size of the "non-flat" regions must be small and x-dependent. Evidence on this point is not 
conclusive. 

Consider now the possibility of strong anisotropy in the Landau function. We see from 
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Eq. |] that because the quasiparticle excitations are concentrated in the nodes, at low T 
the temperature-dependent quantity dL/dT is non-zero only near the zone diagonal. If the 
Landau function f(9\, 62) were anomalously small for either 9\ or 82 near this point then the 
effect of the Landau function on dp s /dT would be greatly reduced in agreement with data. 
The range of angles over which this occurs must be small and f must be x-dependent. 

To summarize: the temperature dependence of the superffuid stiffness strongly suggests 
that some process acting in the high-T c superconductors strongly suppresses the contribution 
of the particles away from the zone diagonal to the current. 
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